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l^ttf"  ABSTRACT 

In  general,    a  reparable  inventory  system  deals  with  two  classes  of 
the  same  type  of  item.      The  first  class  consists  of  those  items  which 
have  been  procured  from  the  manufacturer  and  have  never  been  put  to 
use.      The  second  class  consists  of  those  items  which  have  been  used, 
have  failed,    have  been  repaired,    and  are  ready  to  be  used  again.      In 
addition  to  the  two  classes  of  items,    the  system  also  contains  two  dis- 
tinct inventories.      These  are  the  ready-for-issue  and  the  non-ready-for- 
issue  inventories.      The  ready-for-issue  inventory  contains  both  classes 
of  items  in  their  usable  state,    and  the  non-ready-for-issue  inventory 
contains  only  the  latter  in  its  failed  but  reparable  state.      This  paper 
develops  a  quasi-probabilistic  inventory  model  for  the  reparable  inventory 
system  based  on  the  premise  that  the  repair  facility  is  the  primary  source 
of  inventory,    with  procurements  being  made  periodically  only  to  supple- 
ment this  primary  source. 
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1.       INTRODUCTION 

For  the  most  part,    the  studies  devoted  to  inventory  systems  have 
been  directed  toward  the  consumable  systems.      Recently,    due  to  Defense 
Department  concern  over  the  economics  of  the  systems,    a  considerable 
amount  of  interest  has  been  focussed  on  studies  of  reparable  systems. 
The  studies  to  date  on  the  reparable  systems  have  been  primarily  along 
the  lines  of  a  deterministic  system   [  1  ,  2]  ,    the  interaction  of  the  entities 
of  the  system   [3]  ,    and  glimpses  of  a  probabilistic  system   [4]  . 

The  reparable  inventory  system  differs  from  that  of  the  classical 
textbook  type  in  two  major  ways.      First,    the  system  contains  two  distinct 
inventories,    one  of  which  has  as  its  members  items  that  are  in  a  ready - 
for-issue    (RFI)    state  and  the  other  containing  items  which  are  in  a 
non-ready-for -issue    (NRFI)    state;     i.  e.  ,    the  first  contains  items  which 
are  usable  and  the  second  contains  items  that  must  be  repaired  before 
they  can  be  put  to  use.     Second,    the  RFI  inventory  is  made  up  of  a 
mixture  of  new  items  and  items  that  have  been  used,    failed,    repaired, 
and  are  ready  to  be  used  again. 

The  standard  approach  to  the  problem  of  how  to  run  an  inventory 
system  has  been  to  consider  one  in  which  there  is  a  sole  source  of  supply. 
In  the  reparable  inventory  system,    this  situation  in  general  does  not  exist. 
Here  the  usable  inventory  may  be  thought  of  as  being  supplied  by  two 
separate  processes  which  differ  considerably  from  each  other,    the  first 
being  the  manufacturer  and  the  second  being  a  repair  facility.      These 
processes  may  differ  in  their  leadtimes  and  in  the  cost  of  producing  a 


usable  item.      In  general,    the  repair  leadtime  is  shorter  than  the  procure- 
ment leadtime,    and  the  cost  of  repair  is  less  than  the  cost  of  manufacturing. 
The  arguments  here  tend  to  favor  the  repair  facility  as  a  primary  source 
of  inventory,    but  the  savings  in  operating  this  facility  must  be  large 
enough  to  warrant  its  existence.      The  criteria  that  an  item  must  meet  to 
mark  it  as  reparable  are  threefold.      First,    the  item  must  be  physically 
capable  of  repair;     second,    the  cost  of  accomplishing  this  repair  must  be 
considerably  less  than  the  purchase  price  of  a  new  item;     and  third,    the 
initial  cost  of  the  item  must  be  large  enough  to  warrant  the  cost  associ- 
ated with  getting  the  carcasses  back  to  the  NRFI.     Assuming  that  a 
repaired  item  is  functionally  identical  to  one  which  has  never  been  used, 
it  would  seem  that  the  repair  facility  is  a  more  desirable  supplier  than 
the  manufacturer  and  that  the  decision-maker  should  select  only  the 
facility  to  provide  his  inventory.      This,    in  fact,    would  be  the  case  if  the 
system  were  one  in  which  there  were  no  losses  and  all  carcasses  were 
returned  and  met  the  criteria  of  being  reparable.      For  example,    if  cost 
of  repair  is  the  criterion,    then  the  cost  of  repairing  some  of  the  items 
might  be  considerably  higher  than  the  manufacturing  cost.      In  this 
example,    the  rational  thing  for  the  decision-maker  to  do  would  be  to 
discard  those  items  that  do  not  meet  the  criterion  and  replace  them  by 
procuring  from  the  first  supplier.      Now  assuming  that  the  above  is  the 
case,    then  an  examination  of  the  system  would  indicate  that  a  certain 
fraction   (r)    of  the  items  demanded  from  the  RFI  inventory  will 
eventually  return  to  the  RFI  inventory. 
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Reference   [  1]    deals  with  a  deterministic  model  based  on  a 
"substitution"   policy.      This  policy  requires  the  repair  facility  to  supply 
all  the  items  demanded  until  the  NRFI  inventory  drops  to  a  level  below 
that  necessary  for  another  repair  batch  induction.      The  next  quantity 
that  arrives  at  the  RFI  inventory  is  that  which  had  previously  been  ordered 
through  procurement.     The  demands  are  then  satisfied  with  this  procured 
quantity,    while  the  NRFI  inventory  builds  up  and  the  cycle  starts  again. 
This  paper  will  be  devoted  to  this  substitution  policy,    treating  demand  as 
a  stochastic  variable. 
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2.       NOTATION  AND  ASSUMPTIONS 

The    "substitution"   policy  deals  with  a  system  that  is  self-sustaining 

to  a  certain  point,    at  which  time  a  procurement  must  arrive  in  order  to 

compensate  for  the  losses     ((1    -  r)%)    to  the  system.      Since  the  repair 

facility  is  the  primary  source  of  usable  inventory,    the  model  is  developed 

to  determine  the  most  feasible  batch  size  and  the  number  of  batches  per 

cycle  that  should  be  supplied  by  this  process.      The  model  is  developed 

with  the  measure  of  effectiveness  being  the  cost  of  operating  the  system 

per  year,    using  the  following  notation: 

Q        -     fixed  procurement  quantity,    a  decision  variable; 

Q        -     fixed  repair  batch  size,    a  decision  variable; 
R 

X     -     demand  rate,    units  per  unit  time,    a  random  variable  with 
mean   X   and  density    f  (x  ,  t)    where   t    is  a  particular  time 
period  under  consideration; 

r     -     recovery  rate,    a  given  constant; 

j        -     procurement  leadtime,    a  given  constant; 

T        -     repair  leadtime,    a  given  constant; 
R 

A        -     fixed  cost  of  placing  a  procurement  order,    per  order,    given; 

A        -     fixed  cost  associated  with  an  induction  at  the  repair  facility, 
per  induction,    given; 

h       -     RFI  holding  cost,    dollars  per  unit  per  year,    given; 

h       -     NRFI  holding  cost,    dollars  per  unit  per  year,    given; 

n     -     number  of  inductions  per  cycle,    a  decision  variable, 
n   =    0,    1,    2,    ...      ;' 

II    -     cost  of  incurring  backorders,    dollars  per  unit,    given; 

T        -     expected  time  it  takes  for    Q       items  to  be  demanded; 
R  R 
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6         -     procurement  reorder  level,    based  on  inventory  position, 
a  decision  variable; 

6         -     repair  induction  reorder  level,    based  on  inventory  position, 
a  decision  variable; 

T     -      system  cycle  time,    a  random  variable,    time  between 

successive  procurement  quantity  arrivals  to  RFI  inventory; 

TC        -     total  cost  of  operating  the  system  per  cycle; 

TC     -     total  cost  of  operating  the  system  per  year. 

The    "real"    world  situation  is  a  most  difficult  thing  to  express  in 
terms  of  mathematical  symbols.     Here,    as  in  most  cases,    assumptions 
have  to  be  made  in  order  to  develop  a  model  that,    although  only  an 
approximation  of  the    "real"   world,    gives  an  indication  of  how  the    "real" 
world  behaves.      The  system  considered  in  this  paper  is  one  which 
operates  on  a  cyclic  basis,    where  a  cycle  is  the  random  amount  of 
time    (T)    between  the  arrivals  into  the  RFI  inventory  of  procurement 
orders.      The  cycle  will  consist,    therefore,    of  the  arrival  of  one  procure- 
ment and   n   repair  inductions.     As  long  as  the  distribution  of  the  stochastic 
demand  remains  the  same  from  one  cycle  to  the  next,    an  analysis  of  one 
cycle  will  describe  how  every  other  cycle  behaves.      The  model  presented 
considers  a  system  whose  inventory  is  made  up  of  one  type  of  item.      The 
items  that  are  repaired  are  considered  equal  in  all  respects  to  those 
procured.     The  model  will  be  a  continuous  review-type  based  on  an 
inventory  position,    (  IP  )  ,    where 

IP    =     inventory  on  hand    +     on  order     -    backorders    , 
and  the  on-hand  inventory  is  a  non-negative  quantity.      The  procurement       ■ 
and  repair  orders  are  placed  when  the  IP  drops  to  certain  levels 
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(the  reorder  levels)  .     A  fixed  procurement  order    (  Q     )    is  placed  the 
first  time  the  IP  falls  to  the  procurement  reorder  level  6       after  an 
elapsed  time  of 

n-nr  Q_ 

since  the  arrival  of  a  procurement  order,    where 


T  T 

R 
n       =     

1  TR 


The    [  Z]    is  the  greatest  integer  in    Z.      The  fixed  repair  orders    (  Q     ) 

R 

are  made  the  first  time  the  IP  falls  to  6       after  the  arrival  of  a  repair 
induction.     In  order  to  simplify  the  problem  somewhat,    it  is  assumed 
that  as  a  demand  is  made  on  RFI  inventory,    a  carcass  is  returned  to 
NRFI  inventory  of  which   r%  are  reparable.      With  this  assumption,    it 
is  possible  to  determine  the  exact  amount  of  NRFI  inventory,    given  that 
we  know  the  amount  of  RFI  inventory.     Since  the  system  to  be  formulated 
is  one  in  which  there  is  one  procurement  per  cycle  and  n   repairs  per 
cycle,    the  random  cycle  time    (T)    can  be  expressed  in  terms  of  the 
stochastic  demand   (X)    and  the  quantity  demanded  per  cycle, 

Qp    +    nQ 

T    =  x •  (1) 

In  the  NRFI  inventory,   there  will  be  exactly   rXT   items  that  enter  the 

inventory  each  cycle  and  there  will  be    nQ       items  leaving.     Hence,    in 

R 

steady  state 

rXT    '=    nQR 


14 


or  using  equation  (1)  and  solving  for    Q       yields 

"(1    "  r)  Q 
Qp    *    ; •  (2) 

It  is  postulated  that,    on  the  average,    at  the  time  of  the  arrival  of 
procurement  orders  and  repair  batches  the  net  inventory  will  be  at  fixed 
positive  buffer  levels    b      and   b     ,    respectively.      The  actual  net  inventory 
at  times  of  order  arrivals  will,    of  course,    fluctuate  and,    hence,    the 
buffer  levels  are  indicative  of  protection  against  stockouts. 

Every  time  an  order  is  placed  for  either  a  procurement  or  a  repair 
induction  there  is  associated  a  cost  assumed  to  be  independent  of  the 
number  of  items  either  repaired  or  procured.      These  costs  are  those 
necessary  to  support  the  personnel  and  equipment  involved  in  placing 
orders.      The  cost    (C)    of  the  items  themselves  is  considered  to  be  inde- 
pendent of  the  quantity  procured.      The  RFI  holding  cost   (h    )    is  assumed 
to  be  equal  to  the  cost  of  an  item   (C)    times  an  inventory  carrying  charge 
(I),    i.e., 

h       -     IC   in  dollars  per  unit  year  . 
Since  the  NRFI  inventory  is  made  up  of  just  carcasses,    the  holding  a^ 

cost   (  h    )    is  considerably  less  than  h    .      It  is  assumed  that   h      is  equal 

to  the  cost  of  an  item   (C)  ,   minus  the  cost  of  repair   (C     )  ,    all  times  the 

R 

inventory  carrying  charge    (  I)  ,    i.  e.  , 

h       =     I  (C  -  C     )  ,    in  dollars  per  unit  year     . 
c  R 

The  shortage  cost,    or  cost  of  incurring  backorders,    is  assumed  to 

be  a  function  of  the  number  of  items  backordered.      The  cost  per 
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backorder   (  II  )    is  an  intangible  type  of  cost  and  in  most  cases  should  be 
assigned  by  the  decision-maker.      The  value  of  II   has  a  direct  effect  on 
the  desirability  of  incurring  backorder s,    i.  e.  ,    large  values  of  n   make 
backorders  extremely  undesirable.      The  value  assigned  to   n    is,    in 
most  cases,    much  greater  than  the  value  associated  with  holding  either 
RFI  or  NRFI  inventory.      This  being  the  case,   the  encountering  of  back- 
orders  will  be  the  exception  rather  than  the  rule.     Having  made  this 
assumption,    an  expected  value  formulation,    using  the  expected  values 
as  parameters,    can  be  used  to  determine  the  cost  of  holding  both  RFI 
and  NRFI  inventory;    but  an  exact  procedure  will  have  to  be  employed  to 
determine  the  expected  number  of  backorders.      Figure  1    gives  an 
indication  as  to  how  the  system  really  behaves;    and  Figures  2,3,    and  4 
depict  the  expected  value  formulation. 
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3.       THE  MODEL 

3.  1    Introduction 

With  the  preceding  assumptions  in  hand,    an  expression  will  be 

developed  for  the  expected  cost  of  operating  the  system  per  cycle.     This 

cost,    divided  by  the  cycle  time   ( T)  ,    will  yield  an  expression  for  the 

expected  cost  of  operating  the  system  per  unit  time  in  terms  of  known 

constants  and  the  decision  variables    Q     ,    Q     ,    6     ,   6     ,    and   n  .      The 

P         R        R        P 

cost  per  cycle  is  given  by: 

TC        =     order  cost  per  cycle    (  ORD     )      + 

RFI  holding  cost  per  cycle    (  RFI  HOL     )     + 
NRFI  holding  cost  per  cycle    (  NRFI  HOL     )     + 
shortage  cost  per  cycle    (  SHG_  ) 

The  order  cost  will  simply  be  the  cost  of  placing  one  procurement  order 
plus  the  cost  of  having   n   inductions  from  repair,    i.  e.  , 

ORDT    =    Ap    4-    nAR      .  (3) 

3.  2   RFI  Holding  Cost 

The  expected  holding  cost  per  cycle  for  RFI  inventory  will  be  given 
by  the  product  of  the  holding  cost  per  unit,   per  unit  time    (  h    )  ,    and  the 
area  under  the  RFI  inventory  curve  during  the  cycle.      The  area  under 
the  curve,    see  Figure  2  ,    will  be  the  area  of  the  rectangle  with  sides    b 
and   T_  ,    the  area  of  the  rectangle  with  sides    b      and  nT     ,    the  area  of 

It  C  R 

the  triangle  with  base    T       and  height   Q     ,    and  the  area  of  n  triangles 

with  base    T       and  height   Q      . 
R  R 
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Q      T               n  Q      T 
RFI  HOLT     =    h,  {bL   Tp    +nb2  TR    +    -^-£    +    ^A}     .  (4) 

The  time    T       is  the  expected  time  it  takes  for    Q       items  to  be  demanded; 
hence, 

T        =    

P 

X 

and  (5) 

QR 
R  X 

Substituting  (5)  into  (4)  yields 

2  2 

h  Q  nQ 

RFI  HOLT    =    —     (h1Qp    +nb2QR    +     -|L    +     -y-  }     •  (6) 


Now  using  equation  (2)  , 


n(l    -  r)  Q 
Q. 


P  r 


(4)  becomes 


hl   QR  n        t 
RFI  HOL^    =    - J   b,   (1   -  r)    +    rb. 

T  rX  l     l 

QR    r     2 


r^    [r2    +    n(l-r)2]]     .  (7) 


The  buffer   b      can  now  be  determined  by  again  referring  to  Figure  2 
and  noting  the  inventory  position  just  before  a  procurement  order  is 
placed  and  what  happens  to  the  system  during  the  procurement  lead- 
time   t      .     The  inventory  position  is  at  the  trigger  level  6      .      During 
the  procurement  leadtime    j      ,     all  the  procurement  orders  previously 
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placed  will  have  arrived  and  all  the  outstanding  orders  for  repair  induc- 
tions will  also  have  arrived.     Some    nK    (K  is  to  be  determined)    repair 
inductions  whose  orders  have  not  yet  been  placed  will  also  arrive.      The 
quantity  that  leaves  the  system  during  this  leadtime  will  be  the  mean 
procurement  leadtime  demand     X  t        =     Z     .     Hence,    the  expected  buffer 
will  be 

b,    =    6p    -    Zp    +    „K  QR        .  (8) 

By  similar  analysis,    the  buffer   b      is 

b2    =    6R     "     ZR     "    KQP  '  <9> 

where    KQ       items  previously  ordered  will  not  arrive  during   t      . 
P  R 

Substituting  (8)  and  (9)  into  (7)  yields 

hl  QR  n       / 
RFI  HOLT    =    ^—     J  (6p    -    Zp)    (1    -  r)    +    (ft       -    ZR )  r 

r  X 

QR  r    2 

+ 


2r 


[r       +     n(l    -  r)      ]     J       .  (10) 


In  general,    K  will  be  a  function  of  the  number  of  cycles  of  random   • 
length    T   that  occur  during  the  leadtimes.      Figure  3  ,    for  example, 
indicates  that  there  are  slightly  less  than  three  cycles  occurring  during 
the  repair  leadtime.     Here  it  is  obvious   K    is  equal  to  one.      If   t       re- 
mains the  same  and    Tn    happens  to  be  one  cycle  shorter,    then  in  this 

R 

example     K   =    2    -    0    =    2  ,    which  implies  that  eight  repair  inductions  are 
ordered  and  arrive  during  the  procurement  leadtime.      For  the  general 
case, 
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where    [Z]    indicates  the  greatest  integer  in    Z  . 

3.  3    NRFI  Holding  Cost 

The  NRFI  holding  cost  will  be  obtained  in  a  manner  similar  to  that 

of  the  RFI  holding  cost.      In  this  case,    the  area  in  question  is  depicted 

in  Figure  4.      This  area  will  be  the  areas  of  the  rectangle  with  sides    b 

n(n  -  1) 


and    T,    the  triangle  with  base    T    and  height   nQ      ,      and 

R  4- 

parallelograms  of  horizontal  length    T       and  height   Q      .      The  NRFI 

R  R 

holding  cost  will  then  be 

NRFI  HOLT    =     h2   {b2f    +     £  QRf     -     "'":"    QR2}      •  (12) 


During  the  repair  leadtime,   the  NRFI  inventory  must  supply  the  RFI 
inventory  with    nQ      C     =      J    items.      In  steady  state,    the  level  of  the 


buffer   b      is  the  lowest  possible  level  of  the  NRFI  inventory,    i.  e.  , 


b3 


=    nQR(-T)    =    rZR     •  (U) 


Substituting  equation  (13)  into  (12)  the  NRFI  holding  cost  per  cycle 

becomes 

nQ     T            n(n-  1)   q' 
NRFI  HOLT    =    h2  {  rZRT    +  — }  .    (14) 

3.  4    Shortage  Cost 

Departing  from  an  expected  value  analysis,    an  expression  for  the 
shortage  cost  per  cycle  is  obtained  by  examining  Figures  1  and  2  . 
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Here  backorders  can  be  incurred  at   (n  +   1)    different  times  during  the 

cycle.      There  is  a  possibility  of  stock-out  just  before  the  procurement 

quantity    Q       arrives  and  just  before  the  arrival  of  each  of  the    n   repair 

batches.      The  expected  number  of  shortages,     E  ,    per  cycle  can  be 

obtained  by  examining  what  takes  place  during  the  leadtimes.      Just  before 

a  repair  order  is  placed,    the  IP  is  at  a  level  6     .     As  discussed 

previously,    at  the  end  of  the  repair  leadtime  all  the  items  previously 

ordered  will  have  arrived  with  the  exception  of   KQ     .      The  number  of 

items   short  will  be  zero  as  long  as  the  leadtime  demand  is  less  than 

6D    -   KQ     ;    hence,    the  expected  number  of  items  short  of  the  end  of 
R  P 

n   repair  leadtimes  is 


00 

n 


(x    -   6R    +    KQ     )    f(x;T     )  dx      . 

•5       -KQ  R  P  R 

R  P 

Using  a  similar  analysis  on  the  procurement  leadtime,   the  expected 

number  of  items  short  of  the  end  of  this  period  is 


(x    -   6       -    nKQ     )    f(x;T     )dx      , 
6p    +     nKQp 

where  there  are    nKQ       items  that  are  ordered  and  come  in  during  the 

R 

leadtime.      The  expected  number  of  items   short  per  cycle  will  be 


00 


E     = 


(x    -   6       -    nKQ     )    f  (x;  t     )  dx 
6p+nKQR 


00 

+     n 


(x   -  6_    +    KQD)    f  (x;  t     )  dx     ,         (15) 

aR-KQp 
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where    x   is  a  dummy  variable  indicating  leadtime  demand  and   K    is 
given  by  equation  (11)  .      The  expected  cost  of  backorders  per  cycle  is 
simply 

SHCT    =    IIE       .  (16) 

3.  5   Total  Cost 

The  total  cost  per  cycle  is  given  by  the  expression 

TC       =    ORD^    +    RFI  HOL^    4-    NRFI  HOLm    +    SHC^      . 
T  T  T  T  T 


Since 


(Qp    +    nQR)     t       nQR 


X  - 

r  X 


the  total  cost  per  year  will  be 

TC    =       **        /ORD^    +    RFI  HOL^    +    NRFI  HOL^    4-    SHC^l    .       (17) 
n  O  ^  T  T  1  1 J 


Substituting  equations   (3),    (10),    (14),    (15),    and  ( 1 6)    into  (17),    the  total 
cost  per  year  becomes 
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TC     -  PnQpn     —    *'    +    hl   {(8P    "    V  (1   "  *>    +    <6R    "    Vr 

R 


Q 


R    r    2 


2r 


"r       +     n(l    -  r)  } 


■»{ 


Q 


+     h.         r  Z       + 


R 


R  2 


f  n     -     r  (n  -   1)]  j. 


{J 


4^    J  (x   -  6       -    nKQ     )    f(x;  t     ) 

nQR       lJ6p+nKQR 


+     n 


(x    -6.     +     ^^^    KQJ 
6R    -  i%^    KQR  R  '  R 


i  (  X  ;     TR  )  <lx    \ 


(18) 


Before  attempting  to  minimize  (18)  with  respect  to  the  decision 
variables,  a  value  for  K  must  be  determined.  From  equation  (11), 
K    is  given  as 


K 


zPr 


ZRr 


[nQR  J    "     [  nQR  J 


which  can  be  rewritten  as 


K    = 


zPr 


(^  -DJ  -.Cjfe  -  D0 


Z„  r 
R 


or 


K 


(zp  -  V 


nQ 


R 


(D1     "     °2> 


(19) 
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where 

0    <.    D.    <     1  ,       for     i    =    1  ,   2 


i 


Substituting   (19)    into   (18)  ,     the  total  cost  per  year  becomes 

TC    =  PnV     R       rX    +    h,   {  (6p    -    Zp)   (1  -  r)    ♦    ( 6R   -    ZR)r 

R 

QR    r    2 


♦  -£[■*' ♦■<i-*>'J} 


[»    -    r(n-l)] 


+   TST  1    \      (x-6p  "r(zp  "  zr)  +  nQR(Di  -D2») 

6R+    r(ZF    -    V    -nQR(Dl    "D2» 


f (x;  t     )  dx 


+      n  V  (x   -  6_    +    (1  -  r)  (Z„    -    Z_  )     - 


(20) 
n(l   -  r) 


R  '  v    P  R'  r 


Qr(Dj   -  D2))   f  (x;TR)dx: 


Referring  to  Figure  5,    a  more  exact  analysis  can  be  made  on  the 
values  of   D      and   D     .     It  is  noted  that    D      is  the  time  between  placing 
of  a  procurement  order  and  the  arrival  of  the  next  procurement,    which 
had  previously  been  ordered.     Noting  the  net  inventory,    it  is  apparent 
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that  during  this  same  amount  of  time  some    n    ,  (n     <_    n) ,  inductions  arrive 
from  the  NRFI  inventory.     Hence,      D      can  be  written  as 

Dli    =     "lTRt    SP  TR 
or 


n    r  epr 

D       =      +      

In  n 


where    8       is  that  fraction  between     n     T        and    (n     +   1  )   T        which  is 
also  included  in    D     .     A  similar  analysis  indicates  that    D      is  given  by 

n    r  e_  r 

^                 1  R 

D„     =     + 


2  n  n 


Hence,    the  fraction   D   can  be  written  as 


D=D       -D       =-(e        -6) 
1  2  n    V    P  R' 


From  Figure  5  ,    it  is  readily  apparent  that 


,      „               b2    +     QR     -     S 
€P  TR    <     


X 


or 


b2     -    bl     +     QR 

V    <     Q 

WR 


Also, 


6R    +    QR    "    6P 


CR< 
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Since    €      and  •      are  both  positive,    then 


e       _    e      < 

P  R 


b2     -    bl     -    6R    +    6P 
QR 


where     |  Z  I      is  the  absolute  value  of   Z  .      But, 


bl     "    b2    ~    6P     -    6R 


Consequently, 


€        -    e       ~    o      . 

P  R 


3.  6    Determination  of  Decision  Variables 

Now,    in  order  to  operate  the  system  economically,    values  of  the 

decision  variables  must  be  obtained  which  minimize  the  total  cost  per 

unit  time.      These  values  will  be  obtained  by  setting  the  first  order 

conditions  equal  to  zero.     This  will  be  done  for  the  decision  variables 

Q      ,    6     ,    and  6D     only;    the  integer    n  will  be  determined  by  a  stepping 
RPR 

process,    and   Q       follows    Q      through  equation   (2)  .      This  stepping 
P  R 

process  will  start  by  setting    n   =    1  ,    for  which  a  value  of  Q     ,    6     ,   6     , 

R        R        P 

and    TC    will  be  computed.      The  value  of  n  will  be  increased  by  one, 

and  the  procedure  repeated.      The  process  will  continue  until  the 

minimum  cost  is  obtained.      The  first  order  conditions  are  obtained  by 

taking  the  partial  derivatives  of  the  total  cost  equation   (20)    with  respect 

to  the  decision  variables    Q     ,    6     ,    and  6_    .      First  the  partial  derivative 

R        P  R 

of  total  cost  per  year  with  respect  to   Q       is 

R 
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a  tc 


<Ap    +     nAR) 


3Q 


R 


n  Q 


R 


r  X    + 


4   { 


2  2  i 

r       +    n  (1   -  r)     V 


|n     -     r  (n  -   1)  J 


nrx 


nQ 


R 


(x    -    6n 


r  (Zp    -    ZR))    f  (x;  Tp)  dx 


6P    +     r(ZF    -    V     -     nQR(°l    "    V 


(21) 


+       n 


r 


(x     -    6R    +     (1    -  r)  (Zp    -    ZR))    f  (x;  T     )  dx 


R 


(l-r)(Zp    -    ZR,    +     JU1^IQr,Di-D2)J 


It  becomes  quickly  apparent  that  equation  (21)  cannot  be  solved  explicitly 

for  either     Q     ,    6     ,    or   6      ;    but  an  expression  can  be  obtained  for    Q 
RPR  R 

in  terms  of   Q     ,    6     ,    and  6      .      Setting  equation   (21)    equal  to  zero,    it 
R        P  R 

may  be  written  as 
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Q 


R 


!       \ 


Zr2  X 
n 


h     (r      +    n  (1 


r)     )     +     h2  (rn 


r     (n  -  1)) 


!>    Vap  +  nAR 


)  { 


+    n 


(x    -    6p    -    r  (Zp    -    ZR))     f  (x;  Tp)  dx 


'    V    +     r(ZF    "    V     -     nQR(°l     "    °2) 


(x     -    6R    +     (1    -  r)     (Zp     -     ZR)  ) 
/6R     -     (l-r)(Zp    -    ZR,     +    iLil^lQR(D1-D2) 


f  (x;  TR)  dx 


(22) 


Next,    the  partial  derivative  of  total  cost  per  unit  time  with  respect  to 
6R    is 


a  tc 


56 


R 


V 


IlrX 


Q 


R 


L 


1    -    F(6R     -    (1    -  r)    (Zp    -    ZR) 


♦iiHLV!i 


D2)STR>] 


(23) 


where     F(x;t)     is  the  cumulative  distribution  function  of  demand  during 
time   t  .      Setting  (22)  equal  to  zero,    it  can  be  written  as 
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hlQR 

nx 


[l     -    F,6R    -    (l-r)(Zp        ZR, 


n(l    -  r) 


QR(D1    "    °2);    TR»]    '  ,24) 


Finally,    the  partial  derivative  of  total  cost  per  unit  time  with  respect 
to  6       is 

&  TC  .       ,,  .  IlrX 


h  6n 


=   hi(1"r)   "   T57   [l   ~  f(6p  +  r(zp  "  ZR) 

R 

nQR(Dl    "    °Z);    TP>]      •  (25) 


Here  again  there  is  no  way  of  expressing  any  of  the  decision  variables 
explicitly  in  terras  of  the  other,    but  an  expression  can  be  obtained  which 
will  yield  a  solution  through  iterations  with  equations    (22)    and   (24)  . 
Setting   (25)    equal  to  zero,    it  may  be  written  as 

h     Q      n(l    -r) 
— =     [l     -    F(6p    +    r(Zp-    Z     , 

n  r  X 

"     nQR(Dl    "    D2);TP>]    •  (26) 


If  the  decision  variables    Q       had  been  used  instead  of  Q      , 
equations    (22)  ,    (24)  ,     and   (26)    would  have  been 
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Q, 


2  (1    -  r)2  nX 


2  2  2 

h     (r      +    n(l    -  r)     )    +    h     (rn    -    r     (n  -   1)  ) 


/     ^  A„   +    nAn 
P  R 


+    n 


(x  -  6p  -  r  (Zp  -  Z     ))    f  (x;  Tp)  dx 


^6R+r(Zp    "    V     "    ^TTQP(D1    "    V 


(27) 


+      n 


(x  -   6p   +    (1    -  r)    (Zp    -    ZR))     f  (x;  Tr)  dx 


/ 


'5       -  (1    -  r)    (Z       -    Z     )    +    Q      (D      -    D    ) 
P        v  '    v     P  R'.  P1     1  Z  / 


and 


h,  r  Q_ 
1         P 


n(l    -  r)nX 


=      [l    -    F(6R    -<l-r>(Zp    -    ZR) 

+     QP(D1    "    D2>JTR>] 


(28) 


and 


hiQP 
nx 


=    [l-F(6p+   r(Zp    -    ZR,    -^-Qp,^   -D2>;Tp)].   (29) 


As  a  sidelight,    it  might  be  of  interest  to  note  that  if    r    =    0  ,    above  equations 
would  become 
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Qp 


2X  -l    1/2 

[  —    (Ap   +    n   (  (x    -    6p)    f(x;  Tp)  dx)  J 

1  6„ 


and 


hiQP 
nx 


■  t1    "   f(6p;tp>] 


which  are  identical  to  the  results  obtained  in  consumable  inventory- 
theory   [5]    using  a  similar  approach. 

3.  7    Procedures  to  Obtain  Solutions 

There  are  several  methods  of  obtaining  solutions  to  the  above 
equations.      One  which  makes  no  assumptions  on  the  fraction    D   is  as 
follows: 

1.  Set    n    =    1   . 

2.  Solve  for  an  initial   Q       using  the  equation 

R 


2ApX 
QR  n  (1    -  r)        V  h 


3.       Solve  for         D   =     (D      -    D    )    using  the    Q       from  step  2  as 

12  R 

follows: 

(Z-Z)r  Zr^  Zr 

D    =     -P  R 


CL&J  -LsfcJ) 


nQR  —       -R  -       -R 

where    [  Z]     is  the  greatest  integer  in    Z  . 

4.       Substitute  this  value  of   D  along  with  the  values  of   Q       and   n 

R 

into  equations    (24)    and   (26)    to  determine   6       and  6      . 

R  P 
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5.  Using  the  determined  values  of  Q     ,    6     ,    6     ,    D,    and   n   in  the 

RPR 

right-hand  side  of  (22)  ,    determine  a  new    Q      . 

R 

6.  Continue  steps    3,    4,    and    5   until  the  three  equations  converge. 
If  they  do  not  converge,    assume  a  value  for    D  and  repeat  steps  3 
and  5   until  they  do.     Since  the  value  of   D   is  very  small,    a  good 
value  to  assume  would  be  zero. 

7.  Compute  the  total  cost  per  unit  time  with  the  resulting  values 
of    Q     ,   6     ,   6     ,   n ,    and  K . 

8.  Increase   n   by  one  and  repeat  steps  2  through  8  until  the  minimum 

cost  is  determined.     Select  those  values  of  Q     ,   6     ,   6     ,    and  n 

R       R       P 

which  give  the  minimum  cost  and  compute    Q      from  equation  (2)  . 

A  second  procedure  would  be  to  set    D   =    0   and  repeat  procedure  one, 
keeping   D  at  zero.      This  procedure  should  be  used  only  in  the  event  that 
the  results  of  procedure  one  prove  to  be  infeasible. 


37 


4.       A  NUMERICAL  EXAMPLE 

Since  this  paper  is  an  extension  of  the  substitution  policy,    the 

following  example  is  the  same  as  that  presented  in  reference  [  1]    . 

The  values  of  the  parameters  were  given  as  follows:     A        =     $750  , 

A       =     $100,    r    =    0.  9  ,    X    =     1,000   units  per  year,     h       =     $200  per 
R  1 

unit  year,     h       =     $20  per  unit  year,      t        =     1.0    years, 

w  IT 

t        =     0.  25   years.      The  value  of  h      was  based  on  a  unit  cost   (C)    of 
R  1 

$1,  000  per  unit,    which  implies  a  carrying  charge    (  I)    of  0.  2  per  year, 

and  the  cost  of  repair  of  $900  per  unit.     Here  a  shortage  cost  per  unit 

(H)    is   set  equal  to  $1,  000  ,    and  the  demand  on  the  RFI  inventory  is 

assumed  to  have  a  normal  distribution  with  mean  and  variance  both 

equal  to   Xt  .      The  first  procedure  set  forth  at  the  end  of  the  previous 

section  was  followed  but  the  results,    listed  in  Table  1  ,     obtained  for 

most  of  the  steps  of  the  integer    n   yielded  values  of  the  fraction   D  that 

were  outside  the  permissible  range  set  forth  in  the  previous   section. 

The  value  of  the  fraction  was  set  equal  to  zero,    and  the  procedure  was 

repeated  keeping   D  at  zero.      The  results,    listed  in  Table  2  ,    of  this 

second  procedure  were  almost  identical  to  those  obtained  in  the  first 

except  that  the  reorder  points  changed  considerably.      With  the  fraction 

outside  its  permissible  range,    the  values  obtained  for  the  reorder  points 

were  impossible;    i.  e.  ,    once  the  IP  got  down  to  the   6       level,    it  would 

R 

be  impossible  for  it  to  get  back  up  to  the   6       level.      For  example,    the 

trigger  levels  obtained  for    n   -    17     were     6^    =    356,   6^    =    465,    and 
66  R  P 

the  repair  quantity  was  36  units.      Once  the  IP  got  down  to    6       =    3  56  , 

R 


38 


the  highest  it  could  ever  be  after  this  would  be     6      +    Q       =    392  .     Since 

R  R 

6      +    Q        is  less  than  6      ,     a  procurement  order  would  never  again  be 
R  R  P 

placed  and  the  net  inventory  (NI)   would  tend  toward  negative  infinity. 

The  results  for     D   =    0,    with  rounding  off,    were:     17  lot-size  inductions 

from  the  repair  facility  of  35  units  each;    a  procurement  quantity  of 

67  units;    an  average  cycle  length  of  0.  67    years;    reorder  points,    based 

on  IP,    of    6       =    363    and    6       =    394  ;    all  at  an  expected  total  cost  per 
R  P 

year  of     TC   =    $18,  831.  67  .     It  must  be  mentioned  that  this  total  cost 
does  not  include  the  cost  of  items  nor  the  cost  of  repairing  items.      The 
results  obtained  are  in  excellent  agreement  with  those  in  the  example 
of  reference   [1]  ,    which  are:     19  inductions  of  30  units  each;    a  procure- 
ment of  63  items;    and  a  cycle   (deterministic)   time  of  0.  63    years. 
An  examination  of  Tables  1  and  2    indicates  that  the  buffers   b      and  b 
are  approximately  the  same  for  both  procedures.     With  the  buffers  the 
same,    it  is  expected  that  the  cost  should  be  the  same;     consequently, 
by  setting  the  fraction   D   equal  to  zero,    the  order  quantities  and  costs 
do  not  change  but  trigger  levels  that  are  feasible,    are  obtained.     Again, 
an  examination  of  the  values  in  Tables  1  and  2   quickly  reveals  that  the 
total  cost  of  operating  the  system  is  quite  insensitive  to  the  order  quan- 
tities and  the  number  of  inductions  in  the  neighborhood  of  the  optimum. 
These  results  are  consistent  with  those  obtained  in  consumable  inventory 
theory   [  5]    in  that  a  plot  of  total  cost  versus  procurement  quantity  is 
quite  flat  in  the  neighborhood  of  the  optimum. 
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Even  as  complicated  as  the  equations  appear,    they  are  quite  easily 
solved  with  the  aid  of  a  computer.      The  previous  example  was  programed 
in  FORTRAN  60  and  run  on  a  CDC  1604  ,    with  a  total  run  time  of  slightly 
less  than  two  minutes.      The  results  of  the  iterations  for     n   =    13    and 
D   =    0     are  listed  in  Table  3  .     As  can  be  seen  from  this  table,    the 
equations  converge  quite  rapidly,    which  was  the  case  for  all  steps  of   n 
from  one  through  25  . 

In  order  to  get  some  indication  as  to  how  the  repair  leadtime  affects 

the  total  operating  cost  of  the  system  per  year,    the  above  example  was 

run  again  on  the  computer  with   t       set  equal  to    0.  5  years  vice    0.  25  years, 

R 

The  results  of  this  run  are  listed  in  Table  4.  As  can  be  seen  by  com- 
paring Tables  2  and  4,  the  total  cost  of  maintaining  the  system  can  be 
reduced  for  this  example  by  almost  $8,  000  per  year  by  decreasing  the 
repair  leadtime  from   0.  5  years  to    0.  25  years. 
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5.       CONCLUSION 

In  addition  to  the  basic  assumptions  listed  previously,    the  model 
was  developed  with  one  particular  type  of  reparable  item  in  mind.      This 
type  of  item  can  be  classified  as  being  essential  with  a  high  demand  rate, 
which  usually  implies  that  backorders  are  highly  undesirable.      Items 
that  are  inexpensive,    in  general,    do  not  meet  the  criteria  of  being 
reparable;     hence,    one  more  characteristic  of  the  type  of  item  in  mind 
is  that  it  is  relatively  expensive.      Items  that  cannot  be  handled  by  the 
model  are  those  which  permit  backorders  to  occur  quite  frequently. 
This  class  consists  of  those  items  which  are  quite  expensive  with  a  low 
demand  rate.      To  handle  this  type  of  item,    a  model  would  have  to  be 
developed  using  an  exact  treatment  of  both  the  holding  and  shortage  costs 
Another  type  of  item  that  could  not  be  handled  with  the  model  would  be 
one  in  which  the  distribution  of  demand  changes  considerably  from  one 
cycle  to  the  next.     An  example  of  this  type  would  be  an  item  which  was 
either  being  phased  in  or  being  phased  out  of  the  inventory  system. 

If  the  assumption  dealing  with  the  return  of  carcasses  to  NRFI  is 
beyond  all  realm  of  possibility,    then  the  repair  leadtime  would  have  to 
be  a  random  variable.      This  random  variable  would  be  a  constant   6 
as  long  as  there  are  enough  items  in  NRFI  inventory  to  accommodate 
a  repair  induction  at  the  time  a  repair  order  is  placed.      If  there  are  not 
enough  items  in  NRFI  inventory  at  the  time  an  induction  order  is  placed, 
then  the  leadtime  would  be  increased  a  random  amount  of  time     t    . 
This  would  be  the  time  necessary  to  accumulate  enough  carcasses  to 
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satisfy  the  lot  size  restriction  of   Q       items.      In  this  case,    the  relation- 

R 

ship  between   Q       and  Q       would  not  exist,    and  the  distribution  of  demand 
during  the  repair  leadtime  would  be  a  function  of  two  random  variables. 
Even  with  the  above  type  of  items  excluded,    there  are  still  a  considerable 
number  of  items  that  are  of  the  category  specified  by  the  assumptions, 
and  the  model  presented  should  result  in  values  of  the  decision  variables 
that  are  fairly  accurate. 
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